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Abstract 

A Shashlyk calorimeter prototype for the KOPIO experiment has been con- 
structed and experimentally tested. The energy resolution of about 4%/ v / £(GeV) 
for 0.5 ~r~ 2.0 GeV/c positrons was obtained. Based on this results as well as on 
the results of special measurements, a Monte-Carlo model of the Shashlyk module 
response was developed. This model, including the effects of shower evolution, light 
collection in scintillator plates, light attenuation in fibers, quantum efficiency of 
the photodetector, thresholds and noises in the readout system is consistent with 
experimental results. Possible improvment of the Shashlyk energy resolution up to 
3%/-^/ -E(GeV), the level required by KOPIO experiment, are discussed. 
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1 Introduction. 



The KOPIO experiment [1] at the Brookhaven National Laboratory (BNL) 
Alternating Gradient Sinchrotron (AGS) is designed to measure the decay rate 
for K° L — > tt°i>u, a "gold plated" CP- violating process [2]. This experiment 
will provide the cleanest determination of the fundamental parameter that 
quantifies the phenomenon of CP violation in the context of the Standard 
Model. A measured decay rate very different from the precise expectation of 
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the Standard Model, or one in conflict with CP violation results from the B 
sector, would be evidence for a new physics process. 



A Photon Calorimeter which will occupy an area 5.3 x 5.3 m 2 is one of the 
keystone elements of the KOPIO detector. Studies of Detector optimization 
led to the following requirements for photon detection, in the energy range of 
100 -J- 500 MeV: 



Energy resolution (3 -=- 3.5)%/ '^E7(GeV) 

Time resolution about 80 ps/^/-E(GeV) 
Calorimeter granularity ~ 10 x 10 cm 2 



Shashlyk based calorimeter meets the specified requirements in an economical 
way. Such a calorimeter is composed of Shashlyk modules, which are lead- 
scintillator sandwiches read out by means of Wave Length Shifting (WLS) 
fibers passing through the holes in scintillator and lead. 



The first Shashlyk calorimeter was designed and manufactured at Institute for 
Nuclear Research (Moscow) [3] in 1991 for the experiment 865 [4] (Search for 
the Lepton Number Violating Decay K + — > 7i + fi + e~) at the BNL AGS. During 
the five-year high intensity run of the experiment, the Shashlyk calorimeter 
was a very stable and reliable detector. Its features, together with its low cost 
and well understood method of construction, make this type of calorimeter a 
good candidate for other experimental projects. Similar calorimeters were built 
later for PHENIX experiment at RHIC (BNL) [5] and for HERA-B experiment 
at DESY [6]. A Shashlyk calorimeters was also studied as a candidate for the 
CMS experiment at LHC (CERN) [7], and one is now under construction for 
LHCb experiment [8]. 



The E865 and other constructed calorimeters were designed to have an en- 
ergy resolution about 8%/ \jE(GeV). Significant improvements in calorime- 
ter module construction can be made it possible to achieve the resolution 
~ 3%/y / -E(GeV) required by experiment KOPIO. 



The purpose of this paper is to study the ways of such an upgrade. Based 
on experience with the E865 module, a new prototype module has been de- 
signed, constructed and tested. As will be shown below this module provides 
an energy resolution about 4%/y^(GeV). While this improved resolution is 
not sufficient to meet the requirements of KOPIO, experimental studies of this 
module provide a reference point for tuning Monte-Carlo simulations. From 
such simulations, recommendations for constructing a Shashlyk calorimeter 
for KOPIO are made. 



2 



Scintillator 



WLS fiber 



Fig. 1. Shashlyk module design. 



2 Design of KOPIO prototype module. 



The design of a prototype module for KOPIO is shown in Fig. 1. Eighteen 
such modules had been produced at TECHNOPLAST (Vladimir, Russia). 



2. 1 Mechanical construction of module 



The module is a sandwich of alternating perforated stamped lead and injection 
molded polystyrene-based scintillator plates. The transverse size of module is 
110 x 110 mm 2 ; the thicknesses of plates are 0.35 mm for lead and 1.5 mm 
for scintillator. Each plate has 144 holes equidistantly arranged as 12 x 12 
matrix, with the spacing between the holes being 9.5 mm. The diameters of 
the holes are 1.5 mm in the lead plates, while the holes in the scintillator have 
a conical shape with diameter ranging from 1.4 to 1.5 mm. 72 WLS fibers are 
inserted into these holes. Each fiber is looped at the front of the module, so 
that both ends of a fiber are viewed by the photo- multiplier tube (PMT). Less 
than 10% of light is lost while passing the loop, which radius is about 3 cm. 
The fiber ends are collected in one bunch, cut and polished, and connected 
to a PMT. To improve the light collection 60 /im perforated white reflecting 
paper (TYVEK) is interleaved between lead and scintillator plates, and edges 
of scintillator plates are aluminized. The complete stack of all plates is held in 
compression by the four 50 /im stainless steel side strips that are pretensioned 
and welded to both front and back special endcaps. 

Parameters of module are summarized in Table 1. 



Table 1 

Parameters of the prototype Shashlyk module for KOPIO. 



Transverse size 


110 x 110 mm 2 


Number of the layers 


240 


Polystyrene scintillator thickness 


1.5 mm 


Lead absorber thickness 


0.35 mm 


Reflective paper (TYVEK) thickness 


2 x 0.06 mm 


Number of holes per layer 


12 x 12 


Holes spacing 


9.5 mm 


Holes diameter in Scintillator/Lead 


1.4/1.5 mm 


WLS fibers per module 


72 x 1.3 m rj 92 m 


Diameter of WLS fiber 


1.0 mm, (1.2 mm) 


Diameter of fiber bundle 


14 mm, (17 mm) 


Effective radiation length Xq 


31.5 mm 


Effective Moliere radius Rm 


54.9 mm 


Effective density 


2.75 g/cm 3 


Active length 


473 mm (15.9X ) 


Total length (without photodetector) 


610 mm 


Total weight 


18.0 kg 



2.2 WLS fibers and PM tubes 



We have used three different types of fibers (KURARAY, 1 mm diameter 
Y11(200)M-DC, and BICRON, 1 mm diameter BCF-99-29A-SC and 1.2 mm 
diameter BCF-92-SC) and different types of PMT's (Russian FEU85 and 
FEU115, and 9903B of Electron Tubes Ltd (9903B)) for the tests. 

The properties of WLS fibers used in our measurements are summarized in 
Table 2. These data were taken from manufacturer's Catalogs. Experimental 
measurements [9] of the absorption spectra of WLS fibers in comparison with 
emission spectra of the scintillators are shown in Fig. 2. 

FEU85 is an eleven stages green-extended one inch PM tube with a good 
quantum efficiency for green light (Q.E. ~ 15-20 %). The maximum anode 
current of this PMT is limited to ~ 5 mA per pulse, therefore it was used 
with a fast preamplifier with programmable gain K = 2 4- 10. 

FEU115 is an inexpensive high gain eleven stages PMT with an S20 photo- 
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Table 2 

Properties of WLS fibers. PS stands for polystyrene, PMMA for polymethylmetacry- 
late, and FP for fluorinated polymer. 



WLS fiber 


BCF-92 


BCF-99-29A 


Y11(200)M-DC 


Manufacturer 


BICRON 


BICRON 


KURARAY 


WLS fluor 


G2, 100 mg/1 


G2, 200 mg/1 


K27, 200 mg/1 


Emission peak, nm 


492 


492 


476 


Absorption peak, nm 


410 


410 


430 


Decay time, ns 


2.7 


2.7 


~7 


Core material, 


PS 


PS 


PS 


refractive index 


1.59 


1.59 


1.59 


Inner cladding, 


PMMA 


PMMA 


PMMA 


refractive index 


1.49 


1.49 


1.49 


Outer cladding, 






FP 


refractive index 






1.42 


Trapping efficiency, % 


3.1 


3.1 


5.4 


Attenuation length, m 


~ 3.5 


~ 4.2 


> 3.5 


(for long fiber) 








Fiber diameter, mm 


1.2 


1.0 


1.0 



cathode response and with typical Q.E. ~ 12% for green light. 



9903B is a ten stages green-extended tube with Rubidium bialkali (RbCs) 
photo-cathode and BeCu dynodes for lower gain (0.2 x 10 6 ), but with extended 
linearity (better than 1% up to 50 mA ) and a high short and long term gain 
stability. 9903B has a typical Q.E. ~ 18% for 500 nm. 



2.3 Light collection 

Electrons and positrons from an electromagnetic shower produce ultraviolet 
and blue light in the scintillator plates which is absorbed and reemitted as 
green light in the WLS fiber. About 10% of the green light is captured and 
transmitted by the fiber. 

Due to the Shashlyk design, only light satisfying the criteria of a total internal 
reflection within the scintillator plate is captured by fibers. Since scintillator 
plates are thin compared to the average distance of light to the fibers, many 
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Fig. 2. Absorption spectra of WLS fibers in comparison with emission spectra of 
scintillators. 

reflections from the surface of scintillator occur. Thus, high reflection efficiency 
is very important for good performance. An internal reflection efficiency about 
93% is achievable for realistic surfaces. It should also be noted that modern 
technologies allow comparable diffusion reflection efficiencies to be realized by 
special chemical treatment of the surfaces [11]. 

Light collection uniformity in the direction transverse to the length of the 
module is also important in achieving good energy resolution. This requires 
good reflectivity at the edges of the scintillator plates since the total internal 
reflection at these surfaces is not necessary condition for the possibility of light 
penetrating to fiber. Thus the scintillator edges were aluminized to increase 
light reflection. 



3 Experimental study of prototype module 



The characteristics of the modules were studied on B2 test beam at AGS 
with 0.5 -i- 2 GeV/c positrons and pions. Measurements were done during the 
Spring and Fall test runs 1998. The prototype of calorimeter (3x3 modules) 
was mounted on a platform which could be moved horizontally and vertically 
with respect to the beam line. Upstream of the calorimeter, a trigger counter 
telescope was installed which consisted of three beam defining scintillators 
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(SI - S3) and a 2Xq veto lead-scintillator counter with a 15 mm diameter hole 
through which the beam passed. A Cerenkov counter was used for identifica- 
tion of positrons. 

The signals from PM tubes were read out using an 80 ns wide gate with a 
LeCroy 1885 FASTBUS ADC providing 12 bits of dynamic range. A separate 
channel, including a PM tube without a Shashlyk module and located inside 
the detector box, was used to check for channel to channel coherent noise 
correlations. It was often found that the coherent noise was very high. A 
possible explanation of this effect is that there was no "clean ground" in the 
test beam area. The data from the runs with increased coherent noise were 
skipped in the off-line analysis. 

Three types of triggers were used during the beam tests for monitoring and 
data taking : a random trigger to monitor pedestal behavior and to check noise 
correlations, and two beam triggers (with and without the veto counter) for 
measurements. 

Each module was individually calibrated using either the minimum ionization 
peak from high momentum pions, or the deposited energy peak from 1 GeV 
positrons. In each calibration measurement about 5000 particles were passed 
through the central region of each module at normal incidence. Using these 
data, the calibration coefficients were estimated with precision better than 
1%. 

The energy resolution was determined only for positrons passing through the 
lxl cm 2 area in the center calorimeter as defined by S3 scintillator. 

The momentum spread of the beam significantly contributed to the apparent 
energy resolution in the modules at all measured energies. For beam momen- 
tums above 1 GeV/c, the spread was about 1%. However, according to a 
GEANT [10] calculation, 0.5 GeV positrons passing through matter (15 m of 
Air and 3.5 cm of scintillation counters) lose about 15 MeV of their energy 
with fluctuations of about a ~ 10 MeV. This energy loss required significant 
corrections to the results of the measurements. 

The intensity of beam was chosen sufficiently low to reduce pile-up and rate 
effects during the spill, but due to the large beam size this requirement was 
only partially met. The beam particle rates for a nonet of modules were (50 -j- 
100)K per spill. Small corrections for pile-up and rate effects were taken into 
account in data analysis. 

The electronic noise term for the nonet of modules was measured during spe- 
cial test when calorimeter was removed from the beam. Its contribution to 
the energy resolution of the nonet is 5 ± 1 MeV for FEU85 PM tubes (with 
preamplifiers) and it is 2 zb 1 MeV for other PM's (without preamplifiers). 



7 




0.5 1.0 1.5 2.0 

Energy, GeV 

Fig. 3. The typical pulse height spectrum for the nonet of Shashlyk modules for 
0.5 -=- 2.0 GeV/c positrons. 



The sum of electronic noise and pile-up effect for the nonet was measured in 
other special beam test with the trigger gate for the ADC shifted by about 
300 ns. The contribution of this effective noise to the energy resolution was 
10.3 ± 1.3 MeV for FEU115 PM tubes. This value is somewhat larger than 
expected from pile- up effect, probably due to contributions of "after pulse 
noise" typical for this type of PM tube. The total equivalent noise for FEU85 
PM tubes with preamplifiers was 6.5 ± 1 MeV. The lowest equivalent noise, 
3.6 ± 1 MeV, was obtained with 9903B. 

Signal distributions for 0.5, 1.0, 1.5 and 2.0 GeV/c positrons are shown in 
Fig. 3. Only modules with signals above threshold ~ 3 MeV were included to 
the sum over nonet. We believe that low energy tails for the low momentum 
(0.5 and 1.0 GeV/c) positrons are due to the interaction of the positrons with 
matter upstream of the calorimeter modules. 

The experimental results for energy resolutions for the calorimeter prototype 
with KURARAY Y11(200)-M-DC fibers are shown in Table 3. Corrections to 
energy resolution were obtained by subtracting the contributions of the beam 
positron momentum spread and the effects of energy loss upstream of the 
calorimeter, studied with GEANT Monte-Carlo, from measured values. 

The corrected dependence of the energy resolution on the positron energy may 
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Table 3 

Results of experimental study of Shashlyk prototype (240 layers of 0.35 mm lead 
and 1.5 mm scintillator) with 1 mm diameter KURARAY Y11(200)M-DC fibers. 
Energy resolution is corrected for beam momentum spread and energy loss upstream 
of the calorimeter. 



Pe 


N 




OA 


Measured 


Corrected 


GeV/c 


events 


ADC cnts 


ADC cnts 


a E /E(%) 


a E /E(%) 


2.00 


3000 


1430 


42 


2.92±0.05 


2.70±0.09 


1.75 


3000 


1241 


39 


3.14±0.05 


2.91±0.09 


1.50 


3000 


1056 


35 


3.30±0.05 


3.06±0.10 


1.25 


3000 


913 


34 


3.68±0.07 


3.43±0.11 


1.00 


3000 


699 


29 


4.20±0.07 


3.91±0.11 


0.75 


3000 


531 


25 


4.77±0.09 


4.45±0.13 


0.50 


2500 


359 


23 


6.29±0.17 


5.56±0.20 



be approximated by the following function 



a E , s ~ (3.8 ±0.1)% (-0.8 ±0.6)% 

= (-0.1 ±0.8)%© ; ©- '— ; 0.5<£<2.0, 

E yj E E 

where E is measured in GeV and © stands for qudratic summation. It should 
be emphasized that this dependence was obtained for the nonet of modules 
with a narrow, lxl cm 2 , beam in the center. 

The measured dependence of the energy resolution on energy after correction 
for the positron momentum spread for all studied fiber/PM tubes combina- 
tions is shown in Fig. 4. The GEANT calculated limit of the energy resolution 
due to the sampling only is shown for comparison. 



4 Monte-Carlo simulation 

Energy resolution of a Shashlyk module depends on variety of factors. We 
include the following: 

• Sampling, i.e. thicknesses of lead and scintillator plates. 

• Longitudinal leakage i.e. fluctuation of energy leakage due to the finite 
length of module. 

• Transverse leakage i.e. fluctuation of energy leakage due to the limited num- 
ber of modules used to reconstruct an electromagnetic shower. 

• Effect of the presence of holes, fibers, and steel strips. 

• Light attenuation in the fiber. 
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Fig. 4. The energy resolution of Shashlyk calorimeter (nonet of modules) versus the 
momentum of the positron beam. The structure of calorimeter module is 0.35 mm 
lead, 1.5 mm scintillator, 240 layers. The measurements were done for four types of 
light collection system (see text). 

• Photostatistics. 

• Uniformity of light collection. 

• Electronic noise. 

4.1 GEANT simulation 

GEANT 3.21 was used to simulate the development of electromagnetic showers 
in the module. GEANT contains many tuning parameters which allows one 
to select between speed and quality of simulation. Cuts on energy of electrons 
(CUTELE) and photons (CUTGAM) are crucial for correct simulation of the 
response of Shashlyk modules. Default values for both parameters are 1 MeV. 
The dependence of the energy resolution on the choice of these cuts is displayed 
in Fig. 5. 

Simulated energy resolution is also strongly dependent details of simulating 
energy loss. Default GEANT parameterization does not include generation of 
delta-rays (DRAY=0) and uses Landau- Vavilov-Gauss fluctuations (LOSS=2) 
for energy loss. A significantly different result for energy resolution is obtained 
(Fig. 6) if one enables delta-ray generation (DRAY=1) with appropriate mod- 
ification to energy loss (LOSS=l). 
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Fig. 5. Energy resolution versus GEANT parameter cuts - CUTELE, CUTGAM. 
GEANT simulation was performed for the nonet of Shashlyk modules. The structure 
of a module is 0.35 mm lead, 1.5 mm scintillator, 240 layers. 
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Fig. 6. GEANT simulation of the Shashlyk calorimeter response (en- 
ergy deposited in scintillator) on 250 MeV photon. Dashed line for 
default cuts CUTELE=CUTGAM=1 MeV, DRAY=0, LOSS=2, dotted line 
for CUTELE=CUTGAM=10 KeV, DRAY=0, LOSS=2, and solid line for 
CUTELE=CUTGAM=DCUTE=10 KeV, DRAY=1, LOSS=l. Calculations have been 
performed for long (6000 layers) module consisting of 0.35 mm of lead and 1.5 mm 
of scintillators. Statistics is the same for all three distributions. 

Our calculations use cuts CUTELE=CUTGAM=DCUTE=10 KeV and generat- 
ing of delta-rays (DRAY=1, L0SS=1) since this choice provides better coinci- 
dence with experimental data. 



The dependence of the calorimeter response (visible energy) and energy res- 
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Fig. 7. Calorimeter response and energy resolution as a function of lead and scintil- 
lator plate thickness and number of layers. It is assumed that calorimeter consists 
of 6000 layers of 0.35 mm lead and 1.5 mm of scintillator unless otherwise noted in 
histograms. 

olution on thickness of lead and scintillator plates, and number of layers, are 
shown in Fig. 7. Calculations have been performed for 250 MeV photons. Un- 
less otherwise noted on the histograms, the calorimeter is modeled to be very 
long (6000 layers), consisting of 0.35 mm lead and 1.5 mm scintillator plates. 

Sampling contribution to energy resolution may be estimated by the expression 
a E /VE = 4.7%r 1/2 s 1/4 , where t and s are thicknesses (mm) of lead and 
scintillator plates, respectively. About 20 radiation lengths are required in 
order to achieve a resolution better than 3%/y / i?(GeV), if only the mechanical 
structure of module is taken into account. 



The contribution to the energy resolution of holes in the calorimeter and 
transverse leakage also has been estimated in the GEANT Monte-Carlo and 
will be discussed below. 



4-2 Light attenuation in fibers 



A main concern about attenuation length in fibers is related to the longitudinal 
fluctuations of electromagnetic showers. The typical value for such fluctuations 
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Fig. 8. The effective attenuation of the light in the fibers of Shashlyk module. 
Experimental data (marks) are fit by the exponential dependence exp (—x/L) (solid 
lines), where x is the distance to the PMT and L is the effective attenuation length. 
Measurements have been performed with a muon beam. 



is about one radiation length, which is 3-4 cm for modules of interest (see 
effective radiation length in Fig. 7). This implies that the effective attenuation 
length in fibers installed into the calorimeter modules must be greater than 
2 4- 3 m to have this contribution to energy resolution be much smaller than 
the sampling contribution. 

We have experimentally measured the light attenuation in three different 
fibers: (i) BCF-92-SC, (ii) BCF-99-29A, and (iii) KURARAY Y11(200)M-DC. 
Measurements were performed using muons transversely penetrating modules. 
The size of the beam spot was lxl cm 2 . Results of the measurements are 
presented in Fig. 8. It should be noted that these measurements include ef- 
fects of the fiber loop and the short distance component of light attenuation 
in fiber, i.e., the effective attenuation length in a Shashlyk module for selected 
PM tubes was measured. 

The attenuation length in KURARAY fibers satisfies our requirements. This 
fiber also provides the best light output. For Monte Carlo calculations an 
attenuation length of 300 cm was used, a value which includes the effect of 
fiber loops. 
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Fig. 9. The light yield of the nonet of modules versus the momentum of the positron 
beam. The module structure is 0.35 mm lead, 1.5 mm scintillator, 240 layers. The 
measurements were done for four types of light collection systems (see text). The 
solid lines show the linearity response of detected signals for different energy depo- 
sition in Shashlyk calorimeter. 

4-3 Photo statistics 



Light output was experimentally studied for several types of fibers and pho- 
tomultiplier tubes. The results are displayed in Fig. 9. The best result was 
obtained for KURARAY fiber Y11(200)M-DC. The light output for a module 
with this fiber was experimentally found to be 3100 photoelectrons per GeV. 
This value is matched well by an estimation based on following assumptions: 
(i) 7500 scintillation photons per MeV of deposited energy (for molding scintil- 
lator used in prototype), (ii) 10% efficiency of light collection from scintillator 
to fiber, (iii) 2 x 5.4% trapping efficiency in fiber, where factor 2 corresponds 
to the two directions of light transmission, (iv) attenuation in fiber ~ 0.85, 
(v) quantum efficiency of 9903B photo-cathode ~ 18% 



4-4 Simulation of light collection 



Uniformity of light collection in the scintillator plates is another important 
ingredient for obtaining good energy resolution. 
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To study the dependence of light collection on optical parameters of scintil- 
lator plates, a special optical model was developed. Charged particles passing 
through a scintillator plate generates a number of randomly directed photons 
proportional to energy deposited. Photons so generated may be absorbed on, 
reflect from, or penetrate through the surfaces of the plate. Three types of 
outer surfaces are considered: (i) "clean", i.e., the photon is specularly re- 
flected if the incident angle is appropriate for total internal reflection, other- 
wise it exits from the scintillator, (ii) "mirror", i.e., the photon is specularly 
reflected independently of incident angle, and (iii) "paint", i.e., the photon 
is reflected diffusely. Photons leaving the scintillator may enter a fiber or be 
reflected from environmental "paper" (diffusion reflection) or "mylar" (mir- 
ror reflection). Photons entering a fiber may be reemitted or may exit the 
fiber, depending on the actual length of the photon track in the fiber. The 
model is customized by such parameters as refractive indices of scintillator 
and fiber, attenuation length in scintillator, probability of photon reemitting 
from a fiber, reflection efficiencies, and by geometry (size of the plate, hole 
diameter, fiber diameter). 

In spite of the simplicity of this model (e.g., it does not take into account 
the dependence of reflection probability on incident angle and polarization 
of photons) it gives a reasonable description of light collection in scintillator 
and allows one to compare different choices of the optics of scintillator plate. 
Results of light collection simulation strongly depends on reflection efficiencies, 
which are usually poorly known since they depend on the quality of surface 
treatment. For this reason such a model should be used for relative predictions 
rather than for absolute. 

Model predictions, after tuning of all parameters, for the dependence of light 
output on position of the light source, on the thickness of scintillator plate, 
and on the fiber diameter are compared with experimental measurements in 
Fig. 10-12. 

The following parameters have been used for describing the light collection of 
a prototype module: total internal reflection efficiency equals 93%, the large 
surface of scintillator plate is covered by paper with a reflection efficiency of 
70%, edges of the scintillator plates are mirrored with a reflection coefficient 
of 70%, absorption length in scintillator equals 100 cm, absorption length for 
scintillator light in fiber is 0.5 mm, 5% of light is lost in each passing through 
the boundary between scintillator and fiber. There is an air gap between fibers 
and scintillator. 140 photons per 1 MeV of deposited energy in scintillator were 
generated to match the number of photoelectrons produced by 1 GeV positron 
(if only losses due to the light collection and absorption in fiber are taken into 
account). 
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20 40 

Distance from center of module, mm 

Fig. 10. The nonuniformity of light collection versus the position of the light source. 
The distance is measured from the center of module. The source is a is 1 cm x 1cm 
muon beam. The stars are experimental data, the curve shows the optical model 
simulation. 

5 Monte Carlo study of prototype module 

To check the validity of our Monte-Carlo model, a GEANT simulation of the 
test beam measurements of the prototype modules with 1 mm diameter KU- 
RARAY Y11(200)-M-DC fibers and 9903B was performed. Simulation was 
carried out both for realistic beam, including positron momentum spread in 
beamline and energy loss upstream calorimeter, and for ideal beam. Compar- 
ison of simulation with experimental measurements is displayed in Table 4. 
Experimental and Monte-Carlo with realistic beam results match each other 
with an accuracy of a few percent. This gives us a confidence that Monte- Carlo 
simulation (with an ideal beam) properly reproduces the actual calorimeter 
resolution. 

To understand the contributions of different effects on the energy resolution 
of the module, a Monte-Carlo simulation of the signals produced by 250 MeV 
photons, has been performed. The resulting energy resolutions are displayed 
in Table 5. Sampling effects were simulated in a very long (6000 layers) mod- 
ule. Additional effects were consequently included. Statistical accuracy of the 
calculation is 1.5%. Some disagreement between calculations and test beam 
measurements may be related to the conditions of measurement (beam in 
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2 4 6 

Plate thickness, mm 

Fig. 11. The collection efficiency of scintillation light from scintillator plate for 
green light in WLS fibers versus the thickness of scintillator plate. The circles are 
for experimental data with optical contact between fibers and scintillator, the stars 
are for experimental data without optical contact, and the curves show the optics 
model simulation. 



Table 4 

Experimental and Monte-Carlo mean amplitudes and energy resolutions for proto- 
type Shashlyk module (240 layers of 1.5 mm scintillator and 0.35 mm lead, 1 mm 
diameter KURARAY Y11(200)-M-DC fibers, 9903B). Amplitudes are normalized 
to the 1 GeV amplitude. Experimental energy resolution includes a contribution 
from beam momentum spread, energy loss upstream the calorimeter, and effective 
noise. Only statistical errors are displayed. 

Beam Measurements Monte-Carlo Monte-Carlo 



momentum (realistic beam) (ideal beam) 



(GeV/c) 


<^> 


<r E /E 


(A) 


cte/E 


(A) 


cte/E 


0.5 


0.514 


6.29 ±0.17 


0.487 


5.77 ±0.13 


0.507 


5.17 ±0.08 


1.0 


1. 


4.20 ± 0.07 


1. 


4.22 ±0.08 


1.017 


3.84 ±0.06 


1.5 


1.509 


3.30 ±0.05 


1.514 


3.46 ± 0.06 


1.526 


3.24 ±0.05 


2.0 


2.043 


2.92 ±0.05 


2.031 


3.06 ± 0.09 


2.037 


2.88 ±0.05 
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Fig. 12. Relative light output for different diameter KURARAY Yll fibers compared 
to a 1.0 mm diameter fiber. Experimental points are the data from [11]. The solid 
line is an optical model calculation. 



Table 5 

Calculated energy resolution of prototype module (240 layers of 1.5 mm scintillator 
and 0.35 mm lead, 1 mm diameter KURARAY Y11(200)-M-DC fibers, 9903B). 
uniformly exposed by 250 MeV photons depending on the level of detail of Monte- 
Carlo simulation 

cr E /y/E{GeV) (%) 



Sampling only 2.69 

+ 240 layers 2.84 

+ Holes and steel strips 3.04 

+ Attenuation in fiber 3.01 

+ Photostatistics 3.65 

+ Nonuniformity of light collection 4.07 

+ 3x3 modules 4.27 



the center of the module) in which the effect of the nonuniformity of light 
collection has been suppressed. 
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6 Improving the energy resolution 



One can see that sampling adds the main contribution to energy resolution. 
However, it is less than the combined contribution of other factors 3.3% 
(adding in quadrature). Among these contributions the most significant are 
photostatistics 2.1% and nonuniformity of light collection 1.8%. To reach the 



proposed resolution of 3%/y i?(GeV) sampling, photostatistics and uniformity 
of light collection must be improved. 

The possibilities of improving the sampling contribution are limited. Decreas- 
ing the thickness of lead plates or increasing of thickness of scintillator plates 
increases the effective radiation length of the module, and as a result reso- 
lution deteriorates due to longitudinal fluctuations of electromagnetic shower 
and due to an increased transverse leakage (increased Moliere radius). Si- 
multaneously decreasing the thickness of lead and scintillator plates will lower 
photostatistics. Decreasing lead and/or scintillator plates may also cause tech- 
nical problems for module production. 

To improve the uniformity of light collection we consider using a chemical 
modification of a scintillator surface (CMSS) [11] on the edges of the scintil- 
lator plates, which produces thin (50-100 /zm) white foam layer with a diffuse 
reflection efficiency of about 93%. This coating also increase total light output 
(photostatistics). Monte-Carlo distributions of light output as a function of 
the distance from the center of modules with and without this coating are 
shown in Fig. 13. 

Alternatively, uniformity of light collection may be improved by the appropri- 
ate varying of the density of fiber location in the scintillator plate. 

In addition to increasing scintillator plate thickness and fiber diameter, light 
output may be also increased by establishing optical contact between fiber 
and scintillator. This may be accomplished by gluing fibers within scintillator 
holes. Direct measurements on a single plate (Fig. 11) confirmed this concept. 
The technology of producing this optical contact is not yet developed for a 
total calorimeter, however, so calculations are made for both glue and air 
contact between fiber and scintillator. The gain in photostatistics with glue is 
expected to be about 2. 

Within the limitations on sampling improvement we consider 3 possible ver- 
sions of calorimeter. Results of calculations for the all versions of modules are 
presented in Table 6. 

One can see that improvements in mechanical and optical construction of 
modules can yield energy resolution of about 3%/ JE(GeV), and that modules 
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Fig. 13. Monte-Carlo dependence of light output on distance from center of module 
for (a) existing prototype of module, (b) suggested module with scintillator plates 
totally covered by CMSS, and (c) for module with scintillator plates covered with 
CMSS only on side edge. Distribution are normalized to 1 at center of plate. 

with thinner lead plates might give better results. There is almost no difference 
between the expected resolution for version 1 and 3. A bigger longitudinal 
leakage fluctuations for version 3 are compensated by better sampling term. 

Resolution is approximately 0. 2%/y^GeV) worse if the shower is measured 
by 3 x 3 modules. Increasing of the number of modules to capture a larger 
fraction of visible energy is limited by a greater noise contribution. 

If fibers are not glued, resolution will be additionally increased by 0.1%/ J E(GeV) 



7 Conclusion 



Modules for a Shashlyk calorimeter with energy resolution about 4%/ y E(GeV) 
have been constructed and experimentally tested. Using these results as well 
as results of special measurements, the detailed Monte-Carlo simulation of 
the Shashlyk modules has been developed. Our model of the Shashlyk module 
includes (i) the GEANT 3 simulation of an electromagnetic shower; (ii) the 
optical model for the simulation of light collection in the scintillator plates; 
(hi) effects of light attenuation in the fibers; (iv) effects of quantum efficiency 
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Table 6 

Energy resolution {<je/ \/ E(GeV)) in a module, uniformly exposed by 250 MeV 
photons, depending on the level of detail of the Monte-Carlo simulation. Statistical 
accuracy of the calculations is 2% 



Module version 


1 


2 


3 


Number of layers 


400 


300 


300 


Lead thickness (mm) 


0.25 


0.35 


0.25 


Scint. thickness (mm) 


1.5 


2.0 


2.0 


Fiber diameter (mm) 


1.2 


1.2 


1.2 


Sampling only 


2.14 


2.50 


1.91 


JTJ.I1J.IjC llUIIlUCI Ul IcLyCIh 


9 19 


9 ^4 


9 1 ^ 


+ Holes and steel strips 


2.34 


2.68 


2.31 


+ Attenuation in fiber 


2.53 


2.88 


2.57 


+ Photostatistics 


2.82 


3.09 


2.75 


+ Nonuniformity of light collection 


2.87 


3.15 


2.84 


+ 3x3 modules 


3.17 


3.42 


3.23 


— gluing of fiber 


3.28 


3.53 


3.16 



of the photodetectors, and electronic noise and thresholds in readout system. 
This model gives an excellent description of the energy resolution for variety 
of experimentally studied Shashlyk modules which is illustrated in Fig. 14. 
A few measurements made with low energy photons [12] are also in perfect 
agreement with the developed model. 

Our Monte-Carlo simulation indicates the possibility of improving this resolu- 
tion to about 3%/y / -E(GeV) for module with lead plate thickness 0.25 mm. In 
with a realistic experimental environment, electronic noise and thresh- 
olds, and limited number of modules in a clump, the realized resolution of 
such a calorimeter is expected to be about (3.2 + 3.3)%/y / i?(GeV). 

The experimental study of the test modules was performed in a 0.5+2.0 GeV/c 
positron beam without precise measurement of the beam momentum. For this 
reason extrapolation of the results of the measurements to the KOPIO energy 
region (~ 250 MeV) is somewhat uncertain. The possible way to eliminate this 
ambiguity is a new test measurement in a low momentum electron beam (200+ 
500 MeV/c) with a controlled spread in beam momentum or a measurement 
with low energy photons [12]. 

For future measurements we are planning to use the improved version of the 
KOPIO Shashlyk modules with expected energy resolution (3.0+3.5)%/w E(GeV). 
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• 60 x (4.0 mm Sc + 1 .40 mm Pb) 

T 120 x (3.0 mm Sc + 0.70 mm Pb) 

■ 240 x (1 .5 mm Sc + 0.35 mm Pb) 

v* 300 x (1 .5 mm Sc + 0.25 mm Pb) 
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Fig. 14. Comparison of the Monte-Carlo calculations (solid lines) of the energy 
resolution with our experimental measurements of Shashlyk modules, (•) for BNL 
E865; (▼) prototype for BNL E923; (■) prototype for KOPIO; (*) other prototype for 
KOPIO (not described in this paper). Sampling characteristic, number of layers and 
scintillator /lead thickness is displayed. Filled and empty marks corresponds to the 
measurements with positron and photon beam, respectively. Dashed line indicates 
the expected energy resolution of the upgraded module, satisfying requirements of 
KOPIO experiment. 

These measurements will include a study of the APD photodetectors and Wave 
Form Digitizers for readout. 

Experimental study of time resolution of Shashlyk module is also supposed to 
be done in this measurements. 
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